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Abstract 

This article gives a simple treatment of the quantum Birkhoff nor- 
mal form for semiclassical pseudo-differential operators with smooth 
coefficients. The normal form is applied to describe the discrete spec- 
trum in a generalised non-degenerate potential well, yielding uniform 
estimates in the energy E. This permits a detailed study of the spec- 
trum in various asymptotic regions of the parameters {E, h), and gives 
improvements and new proofs for many of the results in the field. In 
the completely resonant case we show that the pseudo-difFerential op- 
erator can be reduced to a Toeplitz operator on a reduced symplectic 
orbifold. Using this quantum reduction, new spectral asymptotics con- 
cerning the fine structure of eigenvalue clusters are proved. In the case 
of polynomial differential operators, a combinatorial trace formula is 
obtained. 

1 Introduction 

The Birkhoff normal form, in classical mechanics, is a well known refinement 
of the averaging method : under a suitable canonical transformation, a 
perturbation of a harmonic oscillator H2 can be replaced by its average along 
the classical Hamiltonian flow generated by H2- With the averaging method, 
this remains valid as long as one restricts the dynamics to times bounded by 
^(1/e), where e is the size of the perturbation. Using the Birkhoff normal 
form, this time can be extended to for arbitrary N, provided one 

takes into account higher order terms which are also averaged, but in a 
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more intricate sense. Note that, in this work, we do not try to impose 
special restrictions to the original Hamiltonian that would imply some better 
convergence properties (Gevrey convergence, or even analyticity) . Instead, 
we take any smooth function and perform the Birkhoff normal form in a 
neighbourhood of a non-degenerate minimum. 

In quantum mechanics, it is known since at least 1975 that an analogue 
of the Birkhoff normal form can be applied in a very successful way. At 
the formal level, this is attested by physicists like [l^, |l| . Adding on top of 
this the experience of excellent numerical computations, it has become an 
important tool for molecular physics (see [p^] and more recently [20, [2^). 



On the mathematics side, the Birkhoff normal form for pseudo-differential 
operators near a non-degenerate minimum of the symbol has been used by 
several authors already. In particular the article [p^ j by Sjostrand is very 
interesting with this respect, but only deals with the non-resonant normal 
form, that it, when the harmonic oscillator is of the form 

n 

where the coefficients i^j are linearly independent over the rationals. The 
result is that, when the energy E is of order h'^ with 7 G (0, 1) then, for h 
small enough, the quantum system has the same spectrum as a completely 
integrable Hamiltonian. 

The initial goal of our work here is to extend this to the resonant case. 
But since our methods also give new proofs for Sjostrand's result and some 
improvements, and moreover unify them with the analysis of low- lying eigen- 



values initially discovered by Simon |26| and Helffer-Sjostrand [|l^, it might 



be of interest to present it here in the general case. Moreover, we believe 
that several intermediate statements are of independent value, and involve 
for the main part only standard results of semiclassical analysis (symbolic 
and functional calculus for pseudo-differential operators). In particular we 
are able to compare the initial pseudo-differential operator to a differen- 
tial operator with polynomial coefficients, which is very important for many 
practical purposes, including numerical computations. On the other hand, 
the treatment of the resonant case is very hard to perform within the stan- 
dard pseudo-differential calculus because of unavoidable singularities due 
to the fact that, when the coefficients i^j are integers, the S'^-action gener- 
ated by the time-27r flow of H2 is in general non free : periodic orbits with 
smaller periods appear. This explains why so little results were available in 
that case. Our strategy here is to abandon pseudo-differential operators for 
Toeplitz operators, in the spirit of Boutet de Monvel and Guillemin The 
appropriate theory that can deal with orbifold singularities was developed 
in m and il. 
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Let us briefly describe our spectral result in this case. Assume P = 
— + V{x) is a Schrodinger operator with a smooth potential V on 
X = M" or on an n-dimensional compact manifold X equipped with a 
smooth density. (More generally, P could be any pseudo-differential opera- 
tor in some standard class, which is actually our assumption in this article.) 
Assume V G C^{X) has a global minimum at a point which we shall call 
here the origin, and suppose this minimum is non-degenerate. By a linear, 
unitary change of variable in local coordinates near 0, one can always as- 
sume that y'(0) is diagonal; let (z/^, . . . , z/^) be its eigenvalues, with Uj > 0. 
The rescaling Xj ^ y/^Xj transforms P into a perturbation of the harmonic 
oscillator H2: 



P = H2 + W{x), with H2 



where W{x) is a smooth potential of order ^(jxp) at the origin. 

Now assume that the coefficients are completely resonant: there exist 
a real number Uc > and coprime positive integers pi, . . . ,p„ such that 
Vj = Vcpj- Then the spectrum of H2 consists of the arithmetic progression 
En = huci^-^+N) for iV e N, with multiplicity of order A^""^ as A^ ^ 00. It 
is then expected that, for small energies, the spectrum of P is a perturbation 
of the spectrum of H2, splitting each eigenvalue E^ into a band, or cluster. 




We prove this in a precise way. Actually, we describe in theorem 5.3 the 
size and the internal structure of each cluster, as follows. Let H2 be the 
corresponding classical harmonic oscillator, as in (||). It has a 27r periodic 
Hamiltonian flow (pt- Let k = k{x,(^) be the average of W along this flow. 
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k{x,0 = — Wo^{^,{x,()) 

where w : T*X — > AT is the cotangent projection. Let Sn C T*X be the 
sphere: 

Sn = {{x,0(^T*X, H2{x,0 = En}. 
Theorem [theorem 

1. There exists ho > and C > such that for every h G (0, ho] 



Sp{P) n {-00, Chi) c [j 

EjveSp(#2) 



En 3~' ^ ~3~ 



2. When En ^ Chi, let m{EN,h) = #Sp(P) n 



En-!^,En + !^ 



Then m{EN,h) is precisely the dimension ofkei{H2 — En)- 
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3. Let En + \i{En, K),. . . ,Ej\f + \m{Ei^,h){EN ■, K) he the eigenvalues of 
P in this N-eth hand. Then, uniformly for h < Hq and N such that 
En ^ Chi, 

Xi{EN,h)= inf \k{x,C)\ + {EN)-^^{N-^), (2) 
XmiE^,n){EN,h)= sup \k{x,0\ + iEN)^ffiN-^) (3) 

and for any function g £ C°°(M), 

m(Eisi,h) ,/j.x 

where ^iem Liouville measure of Sn 

Thus we see that the average perturbation k behaves as a principal sym- 
bol for the spectral analysis restricted to each cluster. Several improvements 
of this statement are proved in the article. First, k can actually be replaced 
by the homogeneous term of degree 3 in its Taylor expansion. Secondly, the 
exponent 2/3 in the term Ch?/^ (and in {Em)^^"^ where its inverse appears) 
is due to the fact that in general resonances of order 3 may happen in H2 : 
relations of the form pj = 2pj or pj = pj +pfc- If one rules these out, then the 
exponent 2/3 can be replaced by 1/2, but in general with a modified k (if 
the homogeneous term of degree 3 in the Taylor expansion of the potential 
vanishes, then k keeps the same definition. Otherwise the formula is more 
involved). Finally, the last expansion in the theorem is actually the leading 
order of a full asymptotic expansion in h/E. In particular, in sub-principal 
terms, one can exhibit oscillatory contributions of type C,^ where C, is some 
(finite order) complex root of 1. 

The estimate of the spectral density in the particular case pi = ... = 
pn = 1 was first obtained in the thesis of the second author ||2^ through 
a reduction to Toeplitz operators. Independently Bambusi and Tagliaferro 
conjectured and partially proved the estimates of the smallest and largest 
eigenvalues in each band. Then Bambusi and the first author worked on a 
proof using the quantum Birkhoff normal form of |^] and Toeplitz operators. 
The result was announced in 

On technical side, it might be worth mentioning here that we do not use 
any exotic pseudo-differential calculus in order to deal with formal Taylor 
series. Instead of this, we rely extensively on various scaling properties of 
the harmonic oscillator H2, which seems particularly fit for this purpose. 
This allows us to play all along with [E, h) as two (almost) independent 
small parameters. The results of Sjostrand are thus recovered in the regime 
E = W . Then, when we study the resonant case, these scaling properties of 
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H2 become even more crucial, because the effective semiclassical parameter 
becomes h = h/E (sections ^ and P). 

To conclude this introduction, let us mention that quantum Birkhoff 
normal forms have also become a very important tool in inverse spectral 
problems. Formally, the Birkhoff normal form is a (semi)classical invariant 
from which, generally under analyticity assumptions, one can hope to recover 



the full classical dynamics (see for instance |16] and [30|). This aspect is not 
discussed here. 

Structure of the article. — The Birkhoff normal form is based on a 
simple formal construction, which can be explained directly in a quantum 
setting; that's what we recall in section ^ (theorem ^Tll ). Most of the mate- 
rial in this section is not new; however it is crucial here, and the notation 
introduced there is used throughout the article. As we next show in sec- 
tion ^, the relevance of the formal result to semiclassical operators is due to 
a general theorem allowing to compare the spectrum of pseudo-differential 
operators in the so-called semi-excited regime on the basis of the Taylor 



expansions of the symbols (theorem 3^). Adding standard arguments of 



spectral theory we obtain a general statement of the quantum Birkhoff nor- 
mal form (theorem 3.1l| ). Up to an error of size ^{E°°) + &{h°°), it reduces 



the spectral problem to the analysis of a pseudo-differential operator K com- 
muting with a quantum harmonic oscillator i?2- In section ^ we describe the 
joint spectrum of P and its Birkhoff normalisation proving an important 
estimate relating the formal order of K with its operator norm, when re- 



stricted to eigenspaces of H2 (lemma O). Several applications are given : 



Weyl asymptotics, expansions of the low- lying eigenvalues, and the use of 
polynomial differential operators (thus giving a rigorous justification of the 
spectroscopy computations of [^, ^). Finally the last section ^ is devoted 
to the resonant case, when Vj are integers, up to a common multiple. Then 
the classical flow of H2 is periodic, and it is known that one expects the 
spectrum to exhibit clustering. We describe these clusters of eigenvalues. 
Technically and conceptually, the main result is that the restriction of K to 
eigenspaces of H2 can be identified to a Toeplitz operator on the correspond- 
ing reduced symplectic orbifold (theorem |5.lD . This allows us to introduce 
h/E as a. second semiclassical parameter and yields spectral asymptotics for 
these clusters in terms of the principal symbol of K (theorem |5.3| ) . A more 
precise trace formula involving sub-orbifolds and hence oscillatory terms is 



given in 5.5. We end the article with an amusing combinatorial formula 



expressing a certain sum over integral points of a rational polytope, which 



comes as a direct consequence of our results (theorem 5.7 ) 
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2 The formal BirkhofF normal form 

The Weyl quantisation on R^" = T*M" is based on a particular grading for 
formal symbols in x,(,,h, where the degree in the semiclassical parameter 
counts twice the degree of each other variable Xi or ^j. This grading is 
particularly adapted to the harmonic oscillator and hence to the quantum 
Birkhoff normal form. It also appears naturally in the context of deformation 
quantisation |14]. We mainly follow here the presentation of |29], but other 



authors have used this approach. 
Thus we work with the space 

= C|xi,...,x„,^i,...,^„,/i|, 

and we define the weight of the monomial x"'^^h^ to be \a\ + + 2£. The 
finite dimensional vector space spanned by monomials of weight N shall be 
denoted by ^at. Let be the subspace consisting of formal series whose 
coefficients of weight < N vanish. (^Ar)ArgN is a filtration 

^ = ^0 D ^1 D . . . , fl = {0}, 

N 

and shall be used for all convergences in this section. 

The bracket associated to the Weyl product on defines a Poisson al- 
gebra structure on S": it is the unique bilinear bracket for which h is cen- 
tral, which satisfies the Jacobi identity, the Leibniz identity (with the Weyl 
product), and which is commutative on all generators amongst xi, . . . 
^1, . . . and h, except for the relations 

Vj = l,...,re, [^j,Xj] = ^. 

Notice that this structure is invariant by linear canonical changes of coordi- 
nates. There is a simple formula (Moyal's formula) for the brackets of two 
elements of ^, but we shall not need it in this article. However in the fol- 
lowing sections we will use the fact that if H and P are Weyl-quantisations 
of symbols H and P with formal Taylor series at the origin [H] and [P] in 
(f, then the Taylor series of the symbol of the operator commutator [H, P] 
is precisely the Weyl bracket [[H], [P]]. 

The filtration of has a nice behaviour with respect to the Weyl bracket. 
lfNi + N2^2 then 
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If A E the adjoint operator P [^^P] will be denoted by ad^. We shall 
be interested in the adjoint action of elements of Such elements are 
of the form HHq + H2, where Hq G C and H2 is a quadratic form in (x,^). 
Since h is central, we may restrict here to quadratic forms only. They will be 
called elliptic when the quadratic form is positive. Because of this grading 
we see that when H2 £ ^^2, then h^^adn^ acts as an endomorphism of each 
^N- A fundamental property of the Weyl bracket is that ^ad/^jP is exactly 
the classical Poisson bracket {H2,P}. 

We will say that H2 E is admissible when = ker(ad//2)+i™(ad/i-2). 
A typical example is the harmonic oscillator (see lemma below): 

H2 = Mxl + e^)/2 + • • • + i^nixl + fj/2. 

One can show that all elliptic H2 can be written as harmonic oscillators 
in some canonical coordinates and hence are admissible as well. Indeed, 
eigenvalues of Hamiltonian matrices come by pairs (z/j, — z/j) : this implies 
that an elliptic H2 must have the form of a harmonic oscillator plus some 
nilpotent terms. But no such nilpotent term is allowed to show up because 
the flow of H2 is contained in the hypersurface {H2 = const}, which is 
compact. 

The formal quantum Birkhoff normal form can be expressed as follows. 

Theorem 2.1 Let H2 G &2 be admissible and L £ ff^. Then there exists 
Ag and K £ 0'^ such that 

• e'^'^^^{H2 + L) = H2 + K ; 

• [K,H2]=0 . 

Moreover if H2 and L have real coefficients then A and K can be chosen to 
have real coefficients as well. 

Notice that the sum 

is indeed convergent in S' because ^ad^ sends &n into i?N+i- 
Proof. We construct A (and hence K) by successive approximations with 
respect to the filtration of S. Modulo the equality is trivial. So let ^ 1 
and suppose that for some G ^3 we have 

gi^-lad^^ (F2 + L) = + + • • • + Kn+1 + Rn+2 + 

where Ki £ % and commutes with H2, Rn+2 £ S>n+2- Let A' £ !^n+2'-i 
then a small calculation gives 

e'''~"'"'^N+A' ^H2 + L)=H2 + K3 + ---+ Km+1 + Kn+2 + ^7V+3, 
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where 



Kn+2 = Rn+2 + ^adA'H2 = Rn+2 - ^Jt- ^aduiA'. (4) 

We look for an A' such that [i^7v+2) -^^2] = 0. This is possible because H2 is 
admissible. Now if we assume that H2, L and Kj, j ^ N + 1 are real, then 
Rn+2 is real too. Since ^ad/^-j = {H2, •} is a real endomorphism, we have 

^^ = kerK(adH2)+imK(adHj • 

Hence(Q) can be solved with real coefficients. □ 

Remark 2.2 If we write the theorem modulo h we recover the classical 
Birkhoff normal form for Hamiltonians on M^". Indeed, let p and a be C°° 
functions on M^*^ with Taylor expansion at the origin H2{x,^) + L{x,^,0) 
and A{x, ^, 0) respectively. Then if denotes the Hamiltonian flow of a at 
time 1, we have 

po(f) = H2 + k 

where k has the asymptotic expansion K{x,(^,0). Consequently the Poisson 
bracket of H2 and k is flat at the origin. A 

Remark 2.3 Another way of constructing the quantum Birkhoff normal 
form would be to start from the classical result and build successively in 
increasing powers of h. This was used by several authors and amounts 
to follow a different filtration which, in a sense, is less optimal than ours. 



Nevertheless the result is the same, as for instance in |27]. A 



Remark 2.4 The result presented here is often called the Birkhoff-Gustavson 
normal form in the mathematical physics literature. Gustavson popularised 



the idea of Birkhoff in |17| by providing computer programs performing the 
canonical transformation involved. Moreover, Gustavson added the analy- 
sis of the resonant cases, while in his treatise |5|, Birkhoff only dealt with 
the non-resonant situation. Note that Moser had a similar result before 



Gustavson, in the article [22|. In some sense, the Birkhoff normal form is 



the Hamiltonian version of the Poincare-Dulac method |12|. Actually the 
Poincare-Dulac normal form is even more general since it allows for the 
hypothesis of admissibility to be relaxed. Then H2 has to be split into com- 
muting semisimple and nilpotent parts, and the normalisation is performed 
with respect to the semisimple part. The quantum version is probably much 
more complicated to analyse, but it would be very interesting to do so. A 

In this article we will always assume that H2 is elliptic : in some canonical 
coordinates, one can write 



t2^ 
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In order to understand what kind of formal series K can show up in the 
Birkhoff normal form, it is crucial to study the kernel of ad/^j • The following 
lemma is elementary and standard. 

Lemma 2.5 ad/^j is diagonal on the C\h\-hasis where /3,7 are multi- 
indices in N" and Zj = Xj + i^^j , and 

a.dH^{z^z^) = {13- -i,v)zf^zi (5) 

We also state and prove the following — maybe less standard — result, which 
will be one of the tools in the next sections to obtain a pseudo-differential 
version of the Birkhoff normal form. Let be the resonance module 

7^:={aGZ", (a, i/) = 0} (6) 

and denote by n — /c its rank (/c ^ 1). 

Lemma 2.6 There exists a Hamiltonian T'^ action on M^" = r*]R" such 
that the space of all power series that commute with H2 is exactly the space 
of T'' -invariant power series: 

kerad//2 = (T^". (7) 

Proof. Let us use the obvious notation H2 = ^(i^, + ^^). One can 
decompose H2 into 

H2 = l{£i,x''+e)+---+l{ik,x^+e) (8) 

where £j E M.Z", (£1, . . . , are independent, and each Hamiltonian ^{£j,x'^-\- 
has a periodic flow. 

To show this, consider the orthogonal complement of the resonance mod- 
ule 

= {fie Z", (a, I/) = ^ (q, /x) = Va G Z"}. 

Let {u^, . . . ,u^) be a Z-basis of TZ^ . One can view the Q-module = 
7^ (8) Q as a Q- vector space, endowed with the Q-scalar product induced by 
(•,•). Then {TZ^Q)-^ = TZ-^®Q^ and, by density or some algebraic argument, 
(7^ (8> M)-*" = TZ^ (g) M. Therefore v G TZ^ ® M, so we have k ^ 1 and one can 
decompose 

k 
i=i 

We define £j := XjU^ . Since has integer coefficients, it is clear that the 
Hamiltonian ^{u^,x^ + has a periodic flow. 

Prom lemma 2.5 if /?, 7 are multi-indices in N" and Zj = Xj + i^j, then 

s.dH,iz^z'y) = {-f-f3,u)zPz\ (9) 
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But if 7 - /3 G U, then 

7 - /? e Q)^)^ = (7^^ Q)^, 

Hence 

Vj, (n,-,7-/3) =0. 

In other words z^^z'^ commutes with each term in the decomposition (P). 
Therefore any polynomial in keradj^j commutes with all + ■^^)'s, 

and thus is invariant under the T'^ action they generate. □ 

Remark 2.7 Given a point in M^", its orbit under the flow of H2 is con- 
tained in the T'^ orbit of that point. Actually a small variant of the proof 
shows that the inverses of the primitive periods of each periodic Hamilto- 
nian + possess no resonance relation, and hence the //2-01'bit is 

in fact dense in the T'^-orbit. A 

Corollary 2.8 1. If there is no resonance relation (ie. k = n) then any 
element of S commuting with H2 is of the form 

K = f{xl + il---,xl + en:K), 

for a formal series f in n + 1 variables. 
2. More generally if we let 

r = inf{|Q|; q G Z", a / 0, (a, z^) = 0} G N* U {00} 
then any element of S commuting with H2 is of the form 

K = fr{xl+il---,xl+en,K)+Rr, 

where Rr G i^r cind fr{u; h) is a polynomial in (n; h) of degree at most 
[(r-l)/2]. 



Of course this corollary follows even more obviously from lemma 2.5 alone, 
since the monomials z"z^ with |a| + |/3| < r will commute with H2 only if 
a = (3. Hence they admit the form = Jl(^? + 



3 The semiclassical Birkhoff normal form 

The goal of this section is to show how the formal Birkhoff normal form can 
be transformed into a more usable semiclassical statement involving spectral 
estimates. To make the proof more transparent, it is enlightening to separate 
some statements which are independent of the normal form construction, 
and which we believe have their own interest. 
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In all the article we use the following notation. If P is a self-adjoint 
operator on some Hilbert space, P bounded from below, then the increasing 
sequence of eigenvalues below the essential spectrum is denoted by ^ 
AI' ^ • • • ^ AJ" ^ • • • . If / is a borelian of M, the spectral projector of P 
on / is denoted by Uf . UP is a semiclassical pseudo-differential operator, 
then of course A^ = Aj'(/i) and Uf = Ilf{h) also depend on h. 

3.1 Semi-excited spectrum and Taylor expansions 

Let X be either a compact manifold of dimension n equipped with a smooth 
density or X = M". We will deal with semiclassical pseudo-differential 
operators on X in the usual sense, as follows. Let d and m be real numbers. 
When X = W", let S'^{m) = S'^{m,X) the set of all families («(•; /i))?ie(o,i] 
of functions in C'^{T*X) such that 

^C7,n^(l + |x|2 + |C|')^, (10) 

for some constant Ca > 0, uniformly in h. Then ^"^(m, X) is the set of all 
(unbounded) linear operators A on L'^{X) that are ^-Weyl quantisations of 
symbols a G S'^{m) : 

{Au){x) = {Opl{a)u){x) = -i— / ei^^-y^^^a{^-^,i;h)u{y) |«| . 

The number d in (^) is called the order of the operator. Unless specified, 
it will always be zero here. In case X is a compact manifold with a smooth 
density, '^'^{rn,X) is the set of operators on L?{X) that are a locally finite 
sum P = Pp + R, where for each (3 there is a open set Uf^ C X equipped 
with a chart Up — > through which Pp G ^'^(m,M'^), and R is an integral 
operator whose Schwartz kernel is ff{ff^) in the topology. Thus, if 
X is a compact riemannian manifold, A the corresponding Laplacian, and 
V G C°°{X), the Schrodinger operator P = — -^-A-l-yisa good candidate, 
of order zero. In case X = M", the Schrodinger operator is admissible 
whenever V has at most a polynomial growth. 

Let us denote '^{m,X) = Ud(^z^'^{m, X), ^''^(X) = {Jm&z'^'^{m, X), and 
^{X) = Umez^i'iTT', X) . We shall use in this article the standard prop- 
erties of such pseudo-differential operators. In particular the composition 
sends ^{m,X) x ^{m',X) to ^{m + m',X). Moreover all P G ^(0, X) are 
bounded: L'^{X) L'^iX), uniformly for < ^ 1. 

If P has a real-valued Weyl symbol, then it is a symmetric operator on 
with domain C^{X). If its principal symbol is bounded from below then 
we use the Friedrichs self-adjoint extension, and we will identify P with this 
extension. Actually, if the Weyl symbol is real and p is elliptic at infinity (i.e. 
if there exists m G M and a constant C > such that P G ^{m, X) and its 

"I 2 2 / *^ 2 2 

principal symbol p satisfies |p(x, ^)| ^ Tj(||x|| -|- ||^|| )'^' for + ||^|| ^ 



Va G N", 
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C), then P is essentially self-adjoint (see for instance ||Tl|, proposition 8.5]). 
But we won't use this result here. 

Finally, when P G *(m, X) is self-adjoint and / G C^(IR), then f{P) G 
r\mi {fn' , X)) . See for instance or |^ for details. In this work 

all pseudo-differential operators are assumed to admit a classical asymptotic 
expansion in integer powers of h. 

Theorem 3.1 Let P and Q be two semiclassical pseudo-differential opera- 
tors in "^^{X) such that 

• at zq £ T*X , the principal symbols p and q take their minimal value 
p{zq) = q{zQ) = 0, this minimum is reached only at zq and is non- 
degenerate; 

• there exists E^o > such that {p ^ -Eoo} o-nd {q ^ E^o} o^re compact. 

Suppose that, in some local coordinates near zq, the total symbols of P and 
Q admit the same Taylor expansion at zq. Then there exists Eq > 0, Hq > 
and for each N a constant Cm > such that for all {h, E) G [0, Kq] x [0, Eq] 



Af ^E or \^ ^E 



^ Cn{e'' +h 



Before entering the proof of the theorem, we recall an elementary conse- 
quence of the minimax theorem. 

Lemma 3.2 Let A and B be two self-adjoint operators on ^ , both bounded 
from below. Suppose there exists an interval I = {—oo,E] and C > such 
that nf{Jif) C Dom(^) and 

||(^-B)nf II ^ c. 

Then for all j such that X? E one has 

\f ^ Xf + C. 

Proof. Let Af ^ E and = ^S](^) the eigenspace associated to 

the eigenvalues below or equal to A^. Let cj) G , of norm 1; by hypothesis 
one has 

Hence by Cauchy-Schwarz \ {A(j).,(j)) — {B(f),(j))\ ^ C. Therefore 

{Act),ct)) ^ {B(t),4>)+C. 
Since \f = sup^g^s |j^|j^i(5(/>, (/>), one gets 



sup {A(t), 4>) ^ xf + C. 
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Now since J-^ has dimension j, the minimax formula 



Af = inf sup {A(i), <P) 

imphes that \f ^ \f + C . □ 

When deahng with manifolds X ^ M", we shall need a refinement of the 
lemma, as follows. 

Lemma 3.3 Let A and B he two self-adjoint operators acting respectively 
on the Hilhert spaces M" and , both hounded from below. Suppose there 
exists a bounded operator U : — > M" , an interval I = {—oo,E] and 
constants C > 0, c € {0, 1) such that f/Ilf (Jf ) C Dom(yl) and 

\\{U*AU - B)Uf\\ ^ C 

and 

\\u*uuf - nf II ^ c 

Then for all j such that ^ E one has 



Proof. Using the same notation as in the proof of lemma 3.2, we deduce 
from the first hypothesis that 

{AU(j), U(j)) ^ {Be/), <l)) + C, 

while the second yields 



1 



^ c 



Hence \\U Sf > I - c> id. Therefore 



Moreover U : — > J^' is injective and hence dim(C/.Fj^) = j. We conclude 
as in the proof of lemma |3.2|. □ 



Finally, for the proof of theorem 3.1 it will be very convenient to use 



a generalisation of a well-known microlocalisation result of |27[| for which, 
using the above lemmas, we give a new and simple proof. Recall that we 
say that a pseudo-differential operator P £ ^{X) microlocally vanishes at a 
point z £ T*X when in some local coordinates its full Weyl symbol vanishes 
at z. 
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Lemma 3.4 Let P G ^(X') and Q G ^{X) be self-adjoint semiclassical 
pseudo-differential operators, with principal symbols p and q. Assume there 
exists a bounded operator U : L'^{X) — > L'^{X'), compact subsets D C T*X, 
D' C T*X' , and an interval I = (— oo, £J] (with D, D' and E being indepen- 
dent of h) such that 

1. p~^{I) (respectively q^^{I)) is contained in the interior of D (respec- 
tively D' ); 

2. U*PU — Q and U*U — Id are pseudo- differential operators that mi- 
crolocally vanish in D; 

3. P — UQU* and UU* — Id are pseudo- differential operators that mi- 
crolocally vanish in D' ; 



Then there exists > and a positive sequence {Cn)^^^ (depending on 

.f ^ E (or Af 

N 



E ) such that, for all j and h < such that Xf ^ E (or ^ E), one has: 



Remark 3.5 In most situations U will be a Fourier integral operator. 
When X = X' and P, Q G ^(0, X), the result was proved in proposition 
2.2], using the Kato distance to handle the spectral perturbation. Our proof 
here, using the minimax, looks simpler, but the idea is essentially the same. 
A small additional argument is needed to handle X ^ X'. We give the full 
proof here for the sake of completeness. A 

Proof. First let E' > E such that the hypothesis (1.) still holds when E is 
replaced by E' . 

The hypothesis (1.) ensures that the spectra of P and Q are strictly 
bounded from below by some positive constant Eq, independent of h. More- 
over it is well known that it also implies that the intersections with /' = 
{—oo,E'] of these spectra are discrete. One can check this as follows. 

Let El > E' be independent of h and such that q~^{[EQ, Ei]) is contained 
in the interior of D. Let / G C^(M) such that 

1. / =lon [xf,E'] ; 

2. / = outside of [Eo,Ei] . 



By pseudo-differential functional calculus (see for instance [23, theoreme 
III-ll]), f{Q) is a pseudo-differential operator which belongs to the trace 
class and hence is compact. This entails that the spectral projector of Q 
onto [Eq,E'] is compact as well, thus proving the discreteness of Sp(Q) in 
[Eq,E']. Of course the same argument applies to P. See also [jl^ for more 
details. 
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What's more, the functional calculus also ensures that f{Q) microlocally 
vanishes outside q~^{[EQ, Ei]). By symbolic calculus one has 

\\{Q-u*PU)fm\=o{hn. 

Therefore Dom(Qf{Q)) = Dom{U* PU f (Q)) , in the sense of Friedrichs ex- 
tensions. But for all u G Il^,{J^), u = f{Q)u. Hence f/u S Dom(P) and 

\\{Q - U*PU)u\\ = \\{Q- U*PU)f{Q)u\\ = 0{h°°) \\u\\ . 

This shows that there exists a positive sequence {CN)Nm such that for all 



(Q - U*PU)U'f, 



N 



Similarly, \\{U*U - Id)/(Q)|| = 0{h°^) and hence 



([/*[/ -id)n^ 



for a positive sequence {ciy)Nm- 

Applying now lemma we get, for all j such that ^ E' , the in- 
equality Xf ^ (A^ + CNh^){l + iT^^In )- In particular for h smaU enough 
we always reach A J" ^ E' whenever X'j ^ E. 

Interchanging the roles of P and Q we obtainaswell ||(P- C/(3[/*)n|^|| ^ 
Cjsffi'^ and ||([/C/* — Id)n|^|| ^ cnH^ (with perhaps a modification of Cat 
and cat). Hence a new application of lemma ^ yields A^^ ^ (Aj'+CAr/i^)(l-|- 
'^"^^), uniformly for all j such that Xj' ^ E' . This shows that 

P 



-Cm 



A--A^ 



as soon as AV ^ E. Swapping again the roles of P and Q we obtain the 
final result. □ 

Proof of theorem |3.J| . The result of the theorem will be denoted as the 
property Q, ho,Eo, Cn)- It is easy to see that if V{P, Q, Hq, Eq, Cn) and 
V{Q,R,h'Q,E'Q,C'j^) hold, then V{P,R,h!^,E'^,C'l^) will hold with suitably 
chosen constants hQ,EQ,C'^. 

1. — We use this transitivity property to microlocalise the problem in 
a compact subset of T*X. Let <^ be a pseudo-differential operator that is 
microlocally equal to the identity on a neighbourhood of the compact 

:= {p ^ E^} U{q^ E^} 

and microlocally vanishes outside a compact of T*X. We may assume also 
that its principal symbol tp satisfies ^ ^ 1. Then consider the operator 
P' = ^P + 2£^oo(Id — Its principal symbol p' satisfies p' = p for p ^ E^o 
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and p' > Eoo as soon as p > Eoo- Hence P' satisfies the hypothesis of the 
theorem as well. Using that P' is microlocally equal to P on D, we apply 
lemma |3.4| to P and P' with U = Id and an energy E = Eq < E^ such 
that {p ^ E} is contained in the interior of D. Then, since ^ E^ + , 
we obtain V{P, P' , ho, Eq, Cj\[). Using the same trick for Q, we construct Q' 
with V{Q,Q' , h'^, Eq, C'j^) for new constants h'^, Eq, C'j^. 

Thus by transitivity we are reduced to prove the theorem for P' and Q' . 



2. — We compare now P' and Q' . Notice that R ■= P' - Q' = ^{P - Q) 
microlocally vanishes outside a compact subset of T*X. By hypothesis the 
Weyl symbols of P' and Q' near zq have the same Taylor expansion. Hence 
the symbol of R is flat at zq. By symbolic calculus we can construct a 
pseudo-differential operator Sn such that 



and Sn, as R does, microlocally vanishes outside a compact of T*X. This 
implies that Sn is bounded for ?i ^ 1 by a constant independent of h. 



Hence (11) implies, for all ii^ > and h ^ 1, the following estimate 



RU 



p' 



~E,E] 



^ Dn{E 



N , ^N. 



(12) 



We claim that there is a positive sequence (Cat) such that 



RU 



{-oo,E] 



^ Cn{E^ + fi- 



Indeed let — -Emin be the bottom of the spectrum of P' . If E'min ^ E then 



H^' 



n^'^ and the formula follows from (12) with Cn = D^. If 



ii^rnin > E then 



RU 



(-oo,E] 



RUf 



^mim-'^minj 



^ Dn{E, 



N 



But by Garding's inequality there is a constant C > such that -Emin ^ Ch. 
Hence 



RU 



p' 

i-oo,E] 



^ Dn{C^ + l)h'^ ^ Cn{E 



N 



with Cn = Dn{C^ + 1). Thus the claim is proved. 

Now let Eq > Ohe such that the spectrum of P' is discrete in (— oo, Eq]. 
Lemma 3.2 ensures us that if {h, E) G [0, 1] x [0, Eq] then 



Finally, as in the proof of lemma |3.4| , we may interchange the roles of P' 
and Q' to obtain V{P' ,Q' ,hQ, Eq,C'j^), for some new positive constants Hq, 



E'q and C'^. 



□ 
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3.2 Using the formal Birkhoff normal form 

In this section we consider a pseudo-differential operator P £ ^{X) fulfilling 



the hypothesis of theorem 3.1, transform it into a pseudo-differential oper- 
ator on M", take its Taylor series at zq, apply the formal Birkhoff normal 
form, and finally construct a new pseudo-differential operator Q commuting 



with a harmonic oscillator i?2- Q is compared with P using lemma 3.4 and 



theorem 3.1, hence reducing the spectral study of P to that of an effective 



Hamiltonian on some eigenspace of the harmonic oscillator. 
Thus, first of all, we transfer the spectral problem to R"" 



Lemma 3.6 Let P G "^{X) satisfy the hypothesis of theorem 3jJ_ at a point 
Zq G T*X. Then there exists a pseudo-differential operator Q € ^'(M") 
satisfying the hypothesis of theorem \3. l\ at the origin G R^", some constants 
E > and ho > 0, and a positive sequence {Cn)^^^^ such that, for all j and 
h < ho such that ^ E (or ^ E), one has: 



^ CNh 



N 



Moreover there exist local coordinates near zo in which the full Weyl symbol 
of P is exactly the Weyl symbol of Q. 

Proof. Let (x,^) be canonical coordinates on a neighbourhood $7 of = 
(2^0) Co) coming from local coordinates x on X and let U = U{h) be the 
integral operator with Schwartz kernel 

U{x,y) = I ef^<^-y'^^^{y,x,Od^, 

where x G X, y G R" and ip G C^{W^xQ) with (p{y, x,^) = 1 in a neighbour- 
hood of {yo, zo). Then U : L2(M") ^ L'^{X) is bounded and U*PU G ^-CM") 
and has, when expressed in the coordinates (x, C), the same Weyl symbol as 
P. This follows from the fact that, restricted to O and expressed in these 
coordinates, U is simply a compactly supported pseudo-differential operator 
microlocally equal to the identity near the origin. (From a more geometrical 
viewpoint U is actually a Fourier integral operator associated to the sym- 
plectomorphism r*X — > M^" defined by the canonical coordinates (x,^)). 
Since the principal symbol of U*PU has a local non-degenerate minimum at 
the origin, an easy pseudo-differential partition of unity will modify U*PU 
outside a microlocal neighbourhood D of the origin in such a way that its 
principal symbol will satisfy the global hypothesis of theorem Let Q 
be the modified operator. Then for E > small enough, the hypothesis 
of lemma |3.4| are fulfilled (with D' identified with D thanks to the local 
coordinates in 0). This lemma gives exactly the desired spectral result. □ 
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Remark 3.7 We shall not use the fact that P and Q have the same Weyl 
symbols in some coordinates. We only retain the geometrical fact that their 
principal and sub-principal symbols are symplectomorphic. A 

Using the canonical coordinates of T*M"' we introduce the space S' as in 
section |2|. Using a Borel resummation, one can always quantise an element 
in L £ into a pseudo-differential operator in ^'(M"') whose Weyl symbol 
has a Taylor series giving back the initial series in co. Moreover we can 
arbitrarily extend the pseudo-differential operator to vanish microlocally far 
from zq. With a slight abuse with respect to the standard notation, we shall 
in this section denote by Op^{L) such a pseudo-differential operator; and 
for any Q £ ^'(M"') we denote by cr]y{Q) £ S" the Taylor series at zq of the 
Weyl symbol of Q. 

Let Q £ ^{M"-) satisfy the hypothesis of theorem ^ at zq = G M?"". 
Consider the Taylor series [Q] = awiQ)- Since zq is a non-degenerate mini- 
mum for p, one has 

[Q] =hHo + H2 + L 

where L £ ^(3), H2 £ &2 is elliptic, and Hq G M is the value at the origin of 
the sub-principal symbol of Q. Applying the formal Birkhoff normal form 
of theorem O we obtain [A] and [K] in ^3 such that e*'^"''"^!^] {H2 + [L]) = 



H2 + [K]. Consider the operators A = Op^ {[A]) and K = Op^{[K]) (so 
that dwi^) = and awiK) = [K]). Now e*'^ ^ is a Fourier integral 
operator and by Egorov's theorem e*^ '^Qe~^^ £ ^(M"). 

Lemma 3.8 The Taylor expansion at zq of the Weyl symbol of e^^ ^^Qe~''^ 
is hHo + H2 + [K]. 

Proof. Since A is bounded e*'^ ^^Qe~^^ = exp(adj;j-i^)(5. Expanding 
the exponential in the right-hand side using Taylor's formula with integral 
remainder, one gets 

^ 1 1 



where 



Jo 



Rn= I il-tY'{exp{tadif,-iA)Q)dt. 



By definition of the Lie algebra structure of S", <7H/(X]iLo i\,i^'^ifi^^Ay Q) 



precisely J2f=o jr (adj^-i [A] V [Q] = ^Hq + H2 + [K] + ff^+i ■ Thus we need to 
prove that aw i&difi-i Rn) £ ^n+i, and for this it suffices to show 
that the Weyl symbol oi is bounded near zq, uniformly in h. Indeed, its 
Taylor series would then be in and we would conclude using the fact that 
adj;j-i[^] sends ffj to ^j+i- But by Egorov's theorem, exp(tadj;j-i^)Q is a 
pseudo-differential operator of order zero, uniformly in t S [0, 1]. Integrating 
over t we get that Rn is indeed of order zero. □ 
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Proposition 3.9 For any compact D C M^" containing the origin in its 
interior, there exists a pseudo- differential operator K G ^(M""), microlocally 
vanishing outside D, such that 

. [i72,i^] =0; 

• The Weyl symbols of e^^ ^^Qe"*^ and hHo + H2 + K have the same 
Taylor expansion at the origin; 



HHq + H2 + K satisfies the hypothesis of theorem [g. j| (with X 
and zq = 0). 



Proof. We use lemma 2^, ie. the fact that there exists a T action on 
^{xO ('^ith k^l) such that 



kerad^a = (13) 

Let be a compactly supported Borel resummation of [K]. Let K be the 
average of K under the T'^-action. Then, since [K] G (f^*", aw{K) = [K]. 
Hence aw{K) G S''^ and by Weyl quantisation, K := Op^ {K) commutes 
with H2. (Recall that commutation is preserved by Weyl quantisation here 
because H2 is quadratic.) 

Since [K] E ^3, the last point of the proposition is automatically satisfied 
if one chooses the support of the Borel resummation to be close enough to 
the origin. □ 

Prom the proposition, we deduce: 

Corollary 3.10 The operators Q and hHQ + H2 + K have equivalent spectra 



in the sense of theorem 3.1 



Proof. Apply theorem U to e*^ ^^Qe''^ ^ and KHq + H2 + K. □ 

We are now in position to state and prove the quantum BirkhofF normal 
form, which reduces the spectral analysis of P in the semiclassical regime to 
that of a reduced Hamiltonian, ET, acting on some eigenspace of H2 of finite 
dimension (but growing as K decreases or E increases). 

Theorem 3.11 Let P G ^{X) he a semiclassical self-adjoint pseudo-differential 
operator of order zero such that 

• at Zq G T*X, the principal symbol p takes its minimal value p{zq) = 0, 
this minimum is reached only at zq and is non- degenerate; 

• there exists E^o > such that {p ^ -Boo} is compact. 

Let Hq he the value at zq of the sub-principal symbol of P. Then there 
exists a harmonic oscillator H2 on M^" (ie. an elliptic element of &2 in- the 
terminology of section and for any compact domain D C M^" containing 
the origin in its interior there exists a pseudo-differential operator K E 
\I'(M") of order zero such that 
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. [K,F2] =0; 

• K vanishes microlocally outside of D; 

• aw{K) e 

and for each 77 > there exists Eq > 0, Hq > and for each N a constant 
Cat > such that for all {h,E) £ [0,ho] x [0, £'0], 



Af-A? 



hHn 



where 



Q = Q{{1 + V)E) :={H2 + K)^^^^ 



(-oo,(l + r,)B] 



(L2(Rn))' 



Proof. We first apply lemma which allows us to assume that P G ^'(M") 
and zq = 0. 

In view of corollary 3.1C the theorem already holds with Q((l + 'i])E) 
replaced by H2+K. But for E and D small enough the spectra oiQ{{l+rf)E) 
and H2 + K are exactly the same. To see this, we use the following estimate, 
which is due to the fact that [K] G 1^3: 



3^0 > o,3c > 0, > o,v;i < ho 



Kli 



Hi 

i-oo,E'] 



^ CE'^. 



(14) 



This estimate is a particular case of a more general result proven below 
(lemma |4.2D . 

3 

Let Ec = 4/9C^ (this is where the function E' ^ E' - CE'^ reaches 
its maximal value). Then without modifying (|14|) one can assume that 
D is included in H2^{{—oo,Ec]): for this one can replace K by Kf{H2) 
where / E C^(M) takes values in [0, 1], is equal to 1 near the origin, and 
is supported inside {— 00, Ec]. Then the estimate (|l^ can be improved as 
follows: 



VE' > 0, 



^ r{E'), 



(15) 



where r{E') := mm{CE' ^ jCEc^). For any E' > and any operator Q we 
use the notation 



n 



{-oo,E 



Also let Q{oo) := H2 + K. Using that .^^^ ' is stable by H2 and H2 = 
Q{oo) — K we see that 

(16) 



with E = E + 



KIT 



Q(oo) 



^ r{E). Inverting 



But (jl^) implies 

the function E ^ E — r{E) which by construction is strictly increasing for 
^ ^ we see that <^ holds as soon as E ^ E + CeI + ff{E^), which 
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is satisfied E = (1 + ri)E, provided E is small enough. Then ( p^ ) says 
that the eigenvalues of Q(oo) less than E are the same as the eigenvalues 
less that E of the restriction of Q{oo) to J^?^ . This restriction is precisely 
Q{ii + v)E). □ 



4 The joint spectrum 

In the previous section, we have split P into a harmonic oscillator H2 and a 
commuting perturbation K whose microsupport could be arbitrarily small, 
provided one is only interested in sufficiently small eigenvalues of P. As 
m p^, this can be exploited to describe the so-called semi- excited states, 



whose energies are of order 0'(h'^), 7 G (0, 1). Here we push the analysis one 
step further, by introducing a second semiclassical parameter associated 
to the high energy of the harmonic oscillator H2 and which turns out to 
govern the study of the perturbation K. In some sense the idea is to have 
a semiclassical regime associated to the reduction of K by the H2 action. 
This view point is made explicit in the next section where we shall assume 
that this action is periodic. For the moment the game is to control K m a 
semiclassical regime with two semiclassical parameters... 
So, let 

^2(x,o = Ef (17) 

defined on M?^, and let H2 be its Weyl quantisation, acting on L^(]R"). 

Let K = K{h) G ^(M") a pseudodifferential operator commuting with 
H2- Let J^E,h be the eigenspace of H2 for the eigenvalue E. Our goal is 
to study the restriction of K to J^E,h, in terms of the parameters E and 
h, both in a neighbourhood of the origin. More precisely, let E'o > and 
restrict the set of admissible E's to eigenvalues of H2 less that Eq, ie. we 
consider the set 

{{E,h) £ {0,Eo) X {0,00)/ 3a gN", E = h{^4 + i"^"))}- 

Instead of working with {E, h) we shall use the more convenient scaling 
(e, h) defined by 



HE%lifE .€(0,.„U-'6M + (.,N"), (18) 

with eo = ^/Eq. Notice that the dimension of -^Efi can be expressed as a 
function of h alone. Let = We identify ^E,h with using the 

unitary transformation 

U,:^E,h^-yeh, /(x) ^ et/(ex), (19) 
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h h 
Figure 1: Allowed values for the parameters {h,E) or (/i,e). 



The extension of 11^ to L^(M"') satisfies, for any symbol p, 

U,Op^{p{.,h))U: = Opj;M.,h)), (20) 
where q{x,^;h) := p{ex, e^; e^h). 

Definition 4.1 Let K = K{K) he a pseudo- differential operator in ^(R"'). 
For any positive integers I, m we shall say that K S ^i^m if 

• its Weyl symbol k is of order ^{h^) in a neighbourhood of {H2 ^ Eq}; 

• the Taylor expansion of k at the origin is in (with the notation of 
section 

If K{h) commutes with H2, we denote by {K{e,h) : Jifh — > =^)(e,h) the 
family of operators defined by 

The main technical result of this section is the following estimate 

Lemma 4.2 Let K be a pseudo-differential operator commuting with H2. 
Suppose K G ^i^m- Then there exists C > 0, eo > 0, such that 

\\K{e,h)\\^Ce^h'^, VeG(0,eo], V/i"! G + (z^, N"). 

To prove this lemma we shall use a specific version of the functional 
calculus for the harmonic oscillator, which can be obtained as a small im- 
provement of the semiclassical functional calculus of . The only difference 
is that we allow any value of Hq, whereas in [11| /iq has to be chosen "small 
enough". This modification is essential for our purposes. 

Lemma 4.3 Let H2 be a harmonic oscillator as in Then for any 

Hq > and for any f G C^(]R), the family of operators {f{H2))h^ho is 
a semiclassical pseudo- differential operator in ^(0); there exists a bounded 
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family (a(-, in S{0) such that f{H2) = OpJ^{a{-,h)). Furthermore, 
for any integer m, 

a{-, h) = ao + hai + h^a2 H h h^rm{--, h) 

where all ai are Cq° with support in the support of f{H2) (and oq = f{H), 
ai = Oj, and h))fi^ho ^-^ ^ bounded family in S{—N), for all N ^ 0. 

Proof. Assume the result is true for H'q and let us prove it for ho = S'^Hq. 
Let Us be the unitary map (Us'^){x) = 5^'^{5x). Since 

UsH2{h)U*s = 6^H2{h/5^) 

and Us is unitary, we have 

Usf{H2{h))Ul = f{5^H2{h/5^)) 

By the known result, f{6^H2{h')) = Opl,{b{.,h')) for h' ^ h'^. So 

f{H2{h)) = Opl{a{.,h)). ioih^ho 

with a{x,£,,h) = b{6^^x,6~^^, 6~'^h). And we deduce the estimates of a 
from those of b. □ 



Proof of lemma 4J. Let k{x,^;h) be the Weyl symbol of K. First 
(as usual) we can assume that k has compact support. Indeed we split 
k = + where k^ has compact support, /c°° vanishes on 

Be,+s ■.= {z-H2{z)i^Eo + 6} 

for some positive 5, and both k^ and k°^ commute with H2 (using a T'^- 
averaging as in the proof of proposition |3.9D . Then we choose / € C^{M) 
with value 1 on [0, E'o] and support inside (— oo,£'o + (5), so that f{H2) 
and K°° := Op^{k°°) have disjoint microsupport. Then, by lemma 
and symbolic calculus, f{H2)K°° = 0{h^) uniformly for ^ ^ L Since 
h = /^e^ K°°{e,h) is of order &{h^e^) for all N, provided /le^ ^ 2Eo/ W\. 
Now consider q{z; e, h) := k^{ez; e^h). Since k^ has compact support and its 
Taylor expansion is in we have 

q= ^i{^^h)''\ezf) =e^ ^ ^(/i" |z|^). 

2a+l3=e 2a+f3=i 

Since k^ = ff(h^) we know that we can actually write 
q = e^h"' ^{h'^'lzf). 

2a'+f3=i-2m 

We could conclude here by restricting to a compact subset in the new 
variable z and using standard semiclassical arguments similar to what we 
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already used at several places. However this would require to restrict the 
validity of the expected estimate to sufficiently small values of h. So we need 
to refine the analysis. 

By hypothesis h is bounded (by 2/|z^|); so there is a constant C > 
(depending of course on £, m, u) such that 

\q{z-e,h)\ ^ Ce^/i"(z)^-2"', 

For all multi- indices a G N^", d'^k!' still has compact support and its 
Taylor expansion is in So we see similarly that for any a, 

In particular {h~™'e~^Op^{q)) G S{£ — 2m), uniformly for all e ^ eq, h ^ 



2/\iy\. Using lemma O with ho = 2/\i>\ and a function / G C^(IR), we 
consider the operator 

Q := h-"'e-'Opj:{q)f{OpJl{H2)). (21) 

By the composition theorem for pseudo-differential operators (Jill, propo- 
sition 7.7]), Q is uniformly in 5(0), and hence bounded. 

It remains to recall that K{e,h) = Op'^{q)\^^^ which says that 

K{e,h) = h^e'f{l)Q, (22) 

and gives the result. □ 
Let us now turn to some applications of our results. 

4.1 Polynomial differential operators 

If /c G i^i then Op'^ik) is a /i-differential operator with polynomial coeffi- 
cients. Using (^0|) we see that 

U,Opl{k)U: =e'Q{h), (23) 

where Q{h) is a /i-differential operator with polynomial coefficients, inde- 
pendent of e, of order at most i in (x, dx)- 

Now suppose K \s a, ?i-pseudo-differential operator commuting with H2, 
with formal Weyl symbol [K]. Fix £ ^ and write 

with /c^*) G Sii and r^+i G ^{l + 1). Applying lemma O, we get 



K{e, h) = Ko{h) + eKi(/i) + • • • + e^Ki{h) + e^+^Ri+i{h), (24) 
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where is uniformly bounded with respect to h and Ki{h) = Opf^{k^'^^) 

(actuahy, with some abuse of notation, we have written Ki{h) for their re- 
strictions to J^h)- 

In other words, in order to study the spectrum of K{e, h), we can restrict 
to the study of polynomial differential operators commuting with H2, up 
to an error of order ^(e^"*"^), uniformly in h. This is very important in 
practice for numerical calculations, since the matrix elements of such -fCj's 
on a basis of eigenvectors of the harmonic oscillator are very easy to compute 
(especially when the i^j's are written in terms of creation an annihilation 
operators or — this amounts to the same — when studied in the Bargmann- 



Fock representation). See |10| for the case of the 1 : 2 resonance. 

We state the result with the original parameters (E, H) for further refer- 
ence. 

Theorem 4.4 Let P he a h-pseudo-dijferential operator satisfying the hy- 
pothesis of theorem 3.11. Then there exists a harmonic oscillator H2 on M^" 



and a formal series 

[K] = A:(3) + k^^) + . . . G ^3 

commuting with H2 such that, for any i £ 'N and any r] > 0, there exists 
Eq > 0, ho > 0, C > such that for all E G (0, £'0], hG (0, ^0] and j G N, 



Af ^ £ or Xf' ^ E 



Xj' — — hH(j 



^ CE^ + ffih'^), 



where 

Qi := {H2 + Op^ik^^^ + fe(^) + • • • + feW))^ 



Proof. Naturally, we first apply theorem 3.11. Therefore the question is 
reduced to the determination of the spectrum of the operator Q({1 + r])E) 
defined in that theorem, up to an error of order ff{E°°) + 0'{h°°). We need 
to compare Q {{1 + rj)E) to . But equation (U) yields 

Q((l + - mi^,.n = e'+'u:Re+i{h)U, = ff{e'+') 

with e = y/E. Hence Q((l + r])E) - Qi = ri)E)'^ = ff{E^). The 

minimax (lemma |3.2| ) gives the result. □ 

4.2 The semi-excited Weyl law 

Let P be a ?i-pseudo-differential operator satisfying the hypothesis of theo- 



rem 



3.11, and let H2 be the corresponding harmonic oscillator, as in (Q 



We define the resonance order of H2 to be 

:= max(3, inf{|a| ; a G Z", a / 0, (a, i/) = 0}) G N U {00}. 
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By convention = oo is is non resonant. Notice that, in contrast to 
corollary we exclude here the value r,y = 2 since we shall always deal 
with perturbations terms that are formally in ^3. This remark pertains 
through the remaining of the article. 

We are interested here in the counting function 



Theorem 4.5 Let p be the principal symbol of P. For any i < r^ (and 

^P(E,h) = -\- [ \dxdi\ + (^{h-^E^-\h + E'^)) 
uniformly for E and h small enough. 

Corollary 4.6 For any C > 0,7 > 0, £ £ [2,rj/) the Weyl counting function 
for "semi-excited states" has the following asymptotics, as h ^ 0; 



[ZTTtl) Jp^E 



In other words if 7 ^ then the remainder is ^(^~"~'"^'-'^'' 2~)) while if 
7 ^ the remainder is 0'{h^^~'^'^^^~^'^). We see that if 7 < then one 
always gets the best remainder /f(h(^^'^)(^~'y)'' 



Proof of theorem i.t. Apply theorem |4.4| with the same I. Then, using 
the notation of that theorem, 

^P{E,h)=,yV'^'{E,h)+r{E,h), (25) 

where r(S,;i) ^ e [E - Ce'-^ ,E + Ce'-^]]. Since ^ < 

we can see using corollary 2.8 that 

where f{u;h) is a polynomial in {u,h) of degree at most with linear 

part equal to (u, v). The eigenvalues of Qg are therefore of the form 

xQ' = f{h{^ + ai),...,h{^ + an)) 

for integers (ai, . . . , a„) G N", and 

J^Q^{E,,E2; h) := G i^s]} 

= #{aGN"; f {h{\ + ai), . . . ,h{\ + an)) G [^1,^2]}. 
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Using a covering of M" by small cubes with sides of length 2h, centred at 
(^(| + ai), . . . , + an)) , Oi G Z", it is then easy to see that 

nn^Q. (E, ,E2;h)= Vol(/-i ( [E^ ^E^])) 
+ ff{\o\{f-^[Ei -h,Ei + h]))) + ^(Vol(/-i([i^2 -n,E2 + h]))) (26) 

For E small enough, using that f{u;h) ~ (u, z^) with > 0, one has the 
bound 

Yol{r\[E -5,E + 6])) = ^((5^"^^). 
Applying this with 6 = h and 6 = CE~^ we get from ( p6|) 

^'3^(£;-cs^,^ + c^^) = ff{h-''E''--^{h + E^)). 



Another application of (|26D, with Ei = 0, E2 = E, combined with (|2 
yields: 

JK^{E; h) = h-^'Yolif ^ + ff{h-''E''-\h + E^)). 
Notice now that, by a simple Fubini argument, 

Vol{/ i^E} = (2^)-" / \dxdC\ , 



where is the full symbol of Q^. We know from lemma 3.6 and proposi- 
tion U that p = + ^i+i + ff{h). Therefore 



{q^ (^E-C'{h + E^)} c{pi^E}c {q^ i^E + C'{h+E^)}, 
so J^,^j, \dxdC\ = Jp^E \d^dC\ + ffiE^-'ih + E'-^)). □ 

4.3 The low-lying eigenvalues 

Let P be a ^-pseudo-differential operator satisfying the hypothesis of theo- 



rem 



3.11, and let H2 be the corresponding harmonic oscillator, as in ([T7|). 
For an energy E of order the Weyl formula (corollary |4.6| ) says that 
the number of eigenvalues of P below E is bounded, independently of h. 
Of course, this can be obtained directly as a consequence of theorem |4.4| : 
because for E = Ch, the dimension of 11^^^ (L^(M'")) is indepen- 

dent of h. We use here the Birkhoff normal form to recover a result of 
Helffer-Sjostrand [|l^, theorem 3.6] concerning the asymptotics of the low- 
lying eigenvalues of Schrodinger operators. 

Theorem 4.7 For any C > and h small enough, the spectrum of P 
in (—00, Ch] consists of a finite number, independent of h, of eigenvalues. 
These eigenvalues admit an asymptotic expansion of the form 

xf{h) ^ hHo + hfio + d ^hffin,, (27) 

m=0 
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where hfiQ is an eigenvalue of H2, o,nd, as usual, Hq is the value of the sub- 
principal symbol of P at the minimum of its principal symbol. The number 
of eigenvalues with a given //q is equal for h small enough to the multiplicity 
of for H2 ■ In particular the smallest eigenvalue of P is simple for h 
small enough. 



Proof. We apply theorem with E = Ch and £ large. Thus, modulo an 
error of size ff{h^^~^^^^'^), the eigenvalues of P less than Ch are equal to the 
eigenvalues of Qi + HHq. Decomposing the Hilbert space on which Qi acts 
according to the eigenspaces of the harmonic oscillator H2, and using the 
unitary equivalence as in ( p3D , we are reduced to the study of eigenvalues of 
matrices of the form 

Q^(e, h) = e^K^ih) + e^K^ih) + • • • + e^K^ih), 

acting on where h~^ = -^ + (0;, i') for some fixed a £ N" and e = Vhh~^. 
By standard perturbation theory for matrices, the spectrum of e~^Q^(e, h) is 
analytic in e for small e. Since e'^K2{h) is unitarily equivalent to H2 we obtain 
the expansion (p7[), along with the statement concerning the multiplicities. 

□ 

Remark 4.8 In case P is a semiclassical Schrodinger operator —^^lS.-\-V{x)., 



This result appeared almost simultaneously in |26, theorem 5.1] and IT 



theorem 3.6]. The techniques of [19| could in principle be easily generalised 
to treat, like we do here, general pseudo-diff'erential operators. While these 
results are now well known, is it interesting to remark that the appearance 
of half-integer powers of h in the asymptotic expansion was not so obvious 



at that time. Actually in [26| these half-integer powers had been forgotten. 

The Birkhoff' normal form we used here makes it very clear as to why 
and when such powers may appear in the asymptotic expansion of the eigen- 
values. In particular the smallest exponent from which such half-integer 
powers can appear is half the resonance order r^,. Indeed, below this or- 
der, corollary |2.8| shows that the operators Kj(h) must be of the form 

' ' ' ' 2"^" ) fo'^ some polynomial fj and hence have even order in 
(a:;,^). Notice that in dimension 1, r^, = 00; hence only integer powers of 



h may show up in that case. As remarked in [19|, the simplest case where 
half-integer powers of h are present is the so-called 1 : 2 resonance: n = 2 
and V = (1,2). The coefficient of h?/"^ is then the average along the flow of 
H2 of the term of order 3 in the Taylor expansion of the symbol. A more 
general statement is given in section ^.3[ A 



5 Toeplitz operators 



As we saw in equation (|T3|), our initial spectral problem is reduced to the 
spectral analysis of a pseudo-differential operator invariant under a ac- 
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tion, for some k = 1, . . . ,n. With this respect, we have two quaUtatively 
extreme situations: k = n and k = 1. In the first case k = n the harmonic 
oscillator H2 has no resonance relation, and the situation is essentially com- 
pletely integrable, at least in the semi-excited regime. This can be seen 



explicitly with theorem 4.4, in view of corollary 2.5. 

On the contrary, in the case k = 1, the harmonic oscillator H2 is com- 
pletely resonant: up to a common multiple, all frequencies Vi are integers. 
The action is precisely the flow of H2. From the perspective of integrabil- 
ity, the spectral analysis looks more involved. However the fact that H2 has 
a periodic flow forces the spectrum to split into regularly spaced clusters. 
The spectral analysis becomes in some sense simpler, since it comes down to 
the study of each individual cluster. From a geometric point of view each of 
these clusters correspond to the spectrum of an operator acting on a reduced 
space. The goal of this section is to develop this idea. 

In principle, the mixed case 1 < k < n could be treated by a combination 
of both techniques. But this still has to be investigated further. 

We use the notation introduced in the beginning of section ^ In partic- 
ular we deal with rescaled parameters e, h and, using the scaling operator 
Ue as in (1£), we are able to reduce the analysis to the case E = 1. The 
new tool we introduce here is to work in the Bargmann space @] • Recall 
that is the space of entire holomorphic functions on C" with flnite norm, 
where the norm comes from the scalar product 

(V,t/'')C" = / i^,i^'){z)fiiz), with i^,^'){z) = ^(z)^e-l^l'/'^ 

where = ZiZi and fi is the Lebesgue measure on = M?"'. Operators 
on L^(M") can be transported on jS^^ via the Bargmann transform which is 
the unitary map U,gg : L^(M"') given by 

UM^){z) = -4-^ / e^-'i-^'-(^'+^'y')^{x)dx, (28) 
(27r/i)'^"/'i J^n 

where z ■ x = Y^ ZiXi, z"^ = z ■ z, x'^ = x ■ x. 

The harmonic oscillator H2{h) := OpJ^{H2) with H2 as in (jl^) becomes 

" / r) \ 
Hfih) := UaH2ih)mj = hJ2^j [zjg^^ + ij . 

In order to deal with symbols of operators in the Bargmann side we simply 
identify the real phase space M^" with C" using Zj = [xj — Hence 

H2 = T.j^j\zj?- 
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5.1 The reduction setting 

In all this section ^ the main assumption is that for all Vijvj G Q. 
Then there exists positive coprime integers pi, . . . , p„, and a constant z^c > 



such that Vj = Ucpj- Hence, according to (18), our rescaled semiclassical 
parameter h is of the form 

h-^ = Y + ^cN, iV G N. 

Let y = {y G C"; H2{y) = 1}; it is a smooth, compact submanifold of 
C". The Hamiltonian flow of H2 defines a locally free action of on Y 

S^xY^Y, {u,y)^u.y:={u^'yi,...,u^"yn), (29) 

where we identify with the complex numbers of modulus 1. 

The quotient Y/S^ is an orbifold M. It is endowed with the Marsden- 
Weinstein symplectic form lj, naturally defined by 7r*u; = where vr is the 
projection Y ^ M and ivy is the restriction to Y of the canonical symplectic 
form of C". Smooth functions on M are by definition S'^-invariant functions 
on Y. 

Let C^(y) C C°^(Y) be the space of equivariant functions in the fol- 
lowing sense: 

i^{u.y)=u^i;{y), yueS\yeY. 

M is naturally endowed with a complex line bundle L — > y whose sec- 
tions are identified with functions of C^{Y) using the pull-back vr*. More 
generally C°°{M, L^) ~ C^(y). Let us endow L with the hermitian struc- 
ture such that, if 'i/'r and V'r are sections of and ip = 7r*'0r, "0' = ^*^r) 
then 

{A,i^'r)Ky)) = V^(y)V^e-''^^l^l'. (30) 

This defines the scalar product {'tpr,4'r)M '■= fMi''l^r,'4'r)fJ'M, where hm is 
the Liouville (or symplectic) measure of M. 

On the other hand let J^jf^ = keT:{Hf {h) — 1) = U^Jiffi where is, as 
in section the eigenspace of H2{h) corresponding to the eigenvalue E = 1. 
It is well known Q that the monomials z° = • • • z^" such that 

{p,a)=N (31) 



form a basis of =^ • But equation ( plD also shows that the restrictions of 
the z^'s to y belong to C'^{Y). Projecting onto M we thus see that J^i^ 
may be identified with a subspace Hjv C C°°{M, L^). Since in general this 
bijection 

Vn : Hat 
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is not unitary g we introduce Un = {VnY^^Wn : J^jf Hn The 
sequence of spaces (H7v)AreN niust be viewed as the quantising Hilbert space 
for the reduced phase space M, and Un as a kind of Fourier integral operator 
allowing to transport equivariant wave functions of the original Bargmann 
space onto the reduced Hilbert space Hjy. 



5.2 Reduction of the Birkhoff normal form 

Let K = K(h) be an ^-pseudo-differential operator on R" commuting with 
H2- Assume K belongs to the class ^i^rm a-s in definition 4.1, and consider 
again the family of rescaled operators K{e^ h) : J^h — > Since the Taylor 
expansion of the Weyl symbol [K] is we see as in ( p^ ) that the family 
e~^K{e, h) is continuous at e = 0, if we set 

e-'K{e,h)^,=^ = K,{h), 

where Ki(h) is the differential operator with polynomial coefficients obtained 
by the leading term of [K], as precisely defined in (p^). 

We wish to consider K{€, h) as a reduced operator on M . For this 
purpose, let K""^{e,h) = U;jgK{e,h)U^ and introduce 

r(e, h) := e-^UNK'^^ie, h)U% : H;v ^ H^- 

Let IItv denote the orthogonal projector on H^v- 

The main result of this section is that T is a semiclassical Toeplitz op- 
erator on M, as stated in the following theorem. 

Theorem 5.1 If K G ^£,m; then there exists a sequence of functions 
{f{-;N))]\f^fq in C°°{M x [0, eo]) admitting an asymptotic expansion 

fix, e; N) = foix, e) + iV^Vi(x, e) + iV^V2(x, e) + • • • 

for the topology, such that 

Tie, h) = N-'^IlNfi-, e; N) + ff{N~^), 

uniformly in e. Moreover, if km is the K-principal symbol of K then 

e^foi-iriz), e) = k^iez), Vz G Y. 

Proof. The proof proceeds by first representing K as a Toeplitz operator 
on C" using the Bargmann transform and the corresponding formula for the 
Toeplitz symbol. Then one has to show that this symbol can be restricted 
to the energy hypersurface H2 = 1 while retaining the asymptotics in h (or 
N~^) and the invariance. 
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However, some technical preliminaries are required before this. As in 
the proof of lemma L2 we can assume that the symbol k = k{x,^;h) is 
compactly supported. Again we introduce the rescaled symbol 

q{x,^;e,h) = k{ex, e^; e^h). 



Refining the proof of lemma |4?^ we can check that q{-,e,h) can be assumed to 
have a support D that does not depend on e and h and with ^ D. Indeed, 
let $ G Cq°(]R^") with support not containing the origin and identically 
equal to 1 on a a neighbourhood of {H2 = 1}. By /i-symbolic calculus 

ii-op^mip{op^{H2)) = ff(,h^), 

provided the support of if is suitably restricted around 1. Hence if we replace 
q by the Weyl symbol of Op^{q)Op^{^), we deduce from equations (21) 
and (22) that K{e, h) is modified by a term of order ff[e^h°°), which proves 
our claim. Moreover one has 

q{x,^;e,h) = h"'e^q{x,^;e,h) 

where q{-,h) is a family of Cq°(M^" x [0,eo]) admitting an asymptotic ex- 
pansion in powers of h for the topology. Hence, for the proof of the 
theorem, we can safely assume m = i = 0. 

We can now turn to the Bargmann side. Let 

Q^ = UaOp^{q{-,e,h))U*^. 

We know that Q"^ can be represented as a Toeplitz operator. Precisely, 
let 11'* be the orthogonal projector of L^(C",e~l^l onto I^h- For any 
bounded function g on C", the Toeplitz operator with contravariant symbol 
g is by definition the operator 

Then we have the following result 

Theorem 5.2 // g is in the symbol class S{0) on C" ~ M?^ then the op- 
erator U^Tgll^ : L^(]R"') — > L^(M") is a pseudo-differential operator whose 
Weyl symbol is 

^(am = ^ / e-''^-'\^\" g{C + z) \dzd-z\ . 
(vr/i)" Jen 

The map I is continuous from S{0) to S{0). Moreover for any g E S{0) and 
all k ^ 1, 

k-l 

I{9)=Y.(iy^'9/f- + h''Rk{9) (32) 
i=o 

where Rk is a continuous map from 5(0) to 5(0). 
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Proof . Assume first that the symbol g is in the Schwartz class. Since 
U^Uag = Id and UagU*^ = U.^, one has U*jTgU,0g = U^gU,0g and its kernel 
is 

K{x,y)= / U^{x,z)g{z)Up;g{z,y)fi{z) 

11, |2- 



The kernel U;^{z,y) is given in (p^), C/J(x,z) = e ^ 1^1 Uag(z,x) so 

with z = {q — ip)/^/2. We recover the Weyl symbol I{g) with the well-known 
formula 



I{g){x,i)= [ e-^''-'^-^K(x+^,x-^)\dv\ 



2n/2 



(27r/i)3"/2 

= T^ I e-''"^^''+''\{^+P,i + <l)\dqdp\ 

Going back to the coordinates z and C, = {x — i^)/ ^/2, we obtain the desired 
formula for 1(g)- Next following the stationary phase method, we prove that 
the map / is continuous 5(0) 5(0) with the asymptotic expansion (|3^). 
Using a density argument, we conclude that the Weyl symbol of U*jTgU,og is 
I{g) for any g in the class 5(0). □ 

Since the series involved in (|3^) is that of the exponential, it is easy 
to inverse formally. Hence let g{-,h) be a family in C^(C" x [0,eo]) with 
support in D that admits the following asymptotic expansion: 

g{; h) = ^{^yA^q/j\ + ^(/i-). (33) 

Then the proposition says that 

where the remainder is in the uniform norm. Notice that, since is an 
affine function of N with positive slope, asymptotics in h are equivalent to 
asymptotics in A^~^. In particular we can neglect the remainder 0'{h°°). 
Therefore, restricting to M'Jf we can express K""^{e,h) = Q'f^sg as 

K^'ie, h) : ^ 3^^, V ^ g{-. e, (34) 

where H'^i^ is the orthogonal projector of L^(C",e~l^l /'^) onto J^if^ ■ 
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It remains to switch to the space H^v and see how to reduce the symbol 
g = g{-; e, h) to the orbifold M. 

The first step is to prove the theorem with (V^)~^ instead of [/yv In 
other words we look for a suitable symbol J{g) G C°°{M) such that 

or, equivalently, inserting (]3^) , 

U-K^'g = v;jIlNJ{g)VN, 

acting on ■ Using that Yi'^'? and 11^ are self-adjoint, this amounts to 
show that for any ^, ^' G 

{g'>P,^')c" = {J{g)VNi^,VNip')M- (35) 

By definition 

(5V,^')C"=/ e-l^l'/'^<7(^)^(^)?R/"(^), (36) 

where fi is Lebesgue's measure on C" = R^". In order to decompose this 
integral, we use a slicing of C" transversal to Y given by the flow of the 
harmonic oscillator H2 at imaginary times. Precisely, the map 

M X y ^ C" \ {0}, (t, y)^z = (e'.y) := (e'P^yi, . . . , e*P-2/„) 

is a diffeomorphism. Now recall that J^f^ is spanned by (a,p) = A^. So 
any element £ satisfies 

'4}{u.z) = u'^ip{z), Vti G C. 

Hence when n = e* one can write 

|V(e*.z)|'e-'^''l^'-H' = |V'(z)|2e-^^^l-l'e-^^(*'-)e-^h'-H' (37) 

with 

ip{t,z) = -2t + i/, J^(e2P'*- l)|zip. 

Since the origin is not in the support of g we shall now use the coordinates 
{t,y) G R X y in order to calculate the integral ( |36[) . The measure /i can be 
decomposed as 

/i = 6{t,y) \dt\ fiyiy), 

where 5{t, y) is smooth and /xy is the Liouville measure of Y . Recall that 
is 5^-invariant and vr^/xy = ^m- Since /x is also 5^-invariant, the function 
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S{t, y) must be 5^-invariant as well. Finally remark that g also, as defined 
in (|33|), inherits the S'^-invariance of q. This entails, together with (|37|), that 

where the function (g) G C°°{Y) is defined by 

J^(g)(y)= / e~''^^''y^git,y)e"'^\'y\"s{t,y)\dt\. 

Since 99 is 5^ invariant, so is (g). Therefore there exists J{g) G C°°{M) 
such that TT*J{g) = {g). Recall the definition ( pO| ) of the scalar product 
of Htv- Since tt^^^y = /^M> we get the desired identity (|35|). 

Asymptotics of (g) (and hence of J{g)) are obtained by the stationary 
phase lemma. Since ip{t, y) has a global non-degenerate minimum at t = 0, 
the expansion is localised on Y, as expected. For instance, at first order, we 
see that if /o G C°°(M) is such that n* fo{y) = g{0, y; e, h) + ff{h) then 

J{g)im) =e{m)fo{m) + ff{h) 

where e G C°°(M) is a positive function. 

To complete the proof of the theorem it remains to replace (V^)~^ by 
Un- But since Un = {V]\fV^)~^Vpf we have 

UNK-'''{e,h)U*N = {VNVN)'^{V;j)-^K'^{e,h)V^\VNVN)'^, 
which means 

UNK'^{e,h)U*^ = {VnV^)"2 {UNJ{g)) {VnV^)^ . 
If we repeat our argument with K{€, h) = Id +0'{h°°) we see that 

is a Toeplitz operator, and its principal symbol is e. We cannot omit here the 
remainder 0'{h°°) because we work with a compactly supported symbol for 
K. However, one can prove that the uniform norms of iV~4 VXr and its inverse 
are ^(1) (cf. proposition 4.22) Hence (V^)"M^(/i°°))V^^ = ff{h°°) 

and thus is a Toeplitz operator (with asymptotically trivial symbol). Hence 
Wn := {Vn)~^^n^ ^ Toeplitz operator. By the symbolic and functional 
calculus of Toeplitz operator we get that 

U^K'^ie, h)U% = VF^^ {UMJig)) W'^ 
is indeed a Toeplitz operator with principal symbol e~2(e/o)e~2 = /q. □ 
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5.3 Spectral asymptotics of eigenvalue clusters 

In this section we apply the previous result to the BirkhofF normal form of 
a pseudo-differential operator P, in order to get spectral asymptotics that 
we express in the original parameters {E, h). So we assume that P satisfies 
the assumptions of theorem 3.11 with 

n 

For the sake of simplicity, we shall also assume that Hq = : the sub- 
principal symbol of P vanishes at zq. Then, according to this theorem, 
the small eigenvalues of P correspond to eigenvalues of H2 + K, where K 
commutes with H2- 

Let r be equal to 3 when there is a resonance relation of the form pj = 2pi 
or pj = pj + pfc, and to 4 otherwise. Then K belongs to and by lemma 
4.2, the norm of the restriction of K to the eigenspace J^E,h of H2 is ^(£^2 ). 

Since the distance between two consecutive eigenvalues of H2 is larger 
than hi'c, we conclude that the bottom of the spectrum of P split into bands 
or clusters in this precise sense: there exists ho > and C > such that for 
every h £ (0, ho] 

Sp{P)n {-00, Chi) c U [e-'^,e+'^] (38) 

2 

Furthermore for any eigenvalue E for H2 smaller than Chr ^ 
#Sp(P)n[i?-^,^ + ^] =miE,h) 

with the multiplicity m{E, h) = dim J^Efi- First result is an estimate of the 
width of the band and of the distribution of the eigenvalues in each band in 
the regime h/E ^ 0. Denote by 

E + Xi{E, h) ^ .... ^ S + \miE,h){E, h) 

the eigenvalues of P in [E - ^,E + ^]. 
Let ko be the principal /i-symbol of K. 

2 

Theorem 5.3 For any h G (0, ho] and any eigenvalue E ^ Chr for H2, we 

have 

\i{E,h)= inf \ko{x)\+ E^^{h/E), 

xe{H2=E} 

\m(E,n){E,h) = sup \ko{x)\+Ei&{h/E) 

xe{H2=E} 
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and for any function g E C°° (ffi) , 

m{E,h) U ( \ 

where is the Liouville measure of {H2 = E} and the 0"s are uniform 
with respect to h and E. 



To read this result, it is interesting to understand the dependence of the 
leading order terms with respect to h and E as well. Recall that by remark 



2.2 the Taylor expansion of is precisely the classical Birkhoff normal form 
of the principal symbol of P. We have = k^'^'^ -\- i^{\x,^Y'^^) with k^'''\x,^) 
a homogeneous polynomial of degree r. Then 



inf I/cqI = £^2 inf 

{H2=E} {Ii2 = l} 



+ &{E 



2 



L 



- ^ E2 ^ im,=n 



l{H2=E} ^ tj2 y J{H2=l} 

If r = 3, k^'^^ is easily computed. Let the Taylor expansion of the princi- 
pal symbol of P begin with H2 + p^^^ + Then k^^^ is the average 
of p^^^ with respect to the Hamiltonian flow of H2. If r = 4 then in general 
the formula is more involved : letting p = H2 + 0'{\x,^\^) we see 
from the Birkhoff construction that k^'^^ is the average 

where a^^^ is the Hamiltonian whose flow performs the first averaging in the 
method, ie. {H2,a^^^} = —p^^\ We do not know a simpler formulation for 
this term, except, of course, when p^^^ = 0. 

Proof of the theorem. First the result holds for the eigenvalues 

\[{E,h)^...^X'^^j,,^^{E,h) 

of the restriction of K to the eigenspace J^E,h- Indeed, by theorem the 
numbers e~''A^(e^, e^/i) are the eigenvalues of a Toeplitz operator which de- 
pends smoothly on e and with semiclassical parameter h. The principal sym- 
bol of this operator is the push-forward to M of the restriction of e~''A;o(e.) 
to {H2 = 1} or equivalently of the restriction of e~''/co(.) to {H2 = e^}. 
The semiclassical estimates of the smallest and largest eigenvalue and of the 
spectral density of a Toeplitz operator in terms of its principal symbol are 
basic results which were extended to the orbifold case in |^] . Going back to 
the original parameter E and h, we obtain the result for the X^{E,h) and 
the are uniform when E and h run over an arbitrary bounded set. 
Then when E ^ Ch~ and h £ (0, Hq] we have for any 

\X,{E,h)-\'i{E,h)\ ^CN{E^ + h^) 
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and we conclude easily. Since ^Vcfi ^ E ^ Chr the remainders 0'{E^) and 
are negligible when is sufficiently large and disappear. □ 



Remark 5.4 In some sense the result is still true for higher energy E. 
When the condition E ^ Ch~ is violated, the bands may overlap, and we 
can not extract the eigenvalues Xi{E,h) from the whole spectrum of P as 
we did in (^8|). However we saw in the proof that the eigenvalues X'^{E,h) 
satisfy the same estimates as the \i{E,h) except that the remainders are 
uniform with respect to bounded energy E. Therefore the theorem remains 
valid in this regime provided Xi{E, ?i)'s are replaced by X'i{E, ^)'s. Then we 



can recover the spectrum of P with theorem 3.11. A 
K = n * E ~ 00' 

E 
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Figure 2: The spectrum (vertical ticks) exhibits clustering around the eigenvalues 
of the harmonic oscillator (big dots) when E ^ Ch^^^ . But even when overlapping 
occurs, one can give a sense to each band if one thinks in terms of the eigenvalues 



X'i{E, h). See remark 5.4 



The estimate of the spectral density in the previous theorem is actually 
the leading order term of a full asymptotic expansion. The description of 
its various pieces is involved because the reduced phase space M is not a 
manifold but an orbifold. First denote by G the set of complex number (" of 
modulus 1 whose fixed point set 

:= {z G Y; C.z = z] 

is not empty. Here the dot is for the Hamiltonian action generated by H2 
as in (|29|) . A straightforward computation leads to 

G = {C e C*; C^' = 1 for some i] 

and 

Y(^ = C^r\Y with = {z E C"; Zi = if C,^' / 1}. (39) 

The S'^-action preserves Yq. Let be the quotient Y^/S^. It is a twisted 
projective space which embeds into M as a symplectic suborbifold. Denote 
by n{C,) its complex dimension. Finally, let m{C,) be the greatest common 
divisor of {p^; ("P' = 1}. Then with the same argument as in the previous 
proof we deduce from theorem 2.3 of |^] the 
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Theorem 5.5 For every function g G C°°(M), 

'"E.(A.(...)/.i)^EQ"'V"f;(|)Vi.o.^((|)°°) 

where the remainder is uniform when h and E runs over (0, Kq) and Sp(ff2)n 

2 

{ — oo, Ch~) respectively, and N is defined by E = + N). Furthermore 
the coefficients Ii{€,C) are C°° function of e and 

where is the Liouville measure of M^. 



5.4 A convex polytope and a trace formula 



As a consequence of the preceding theorems |5.3| and 5.5 we obtain an in- 
teresting formula expressing the asymptotics of a combinatoric sum over 
integral points of a convex polytope, when the polytope undergoes some 
rational scaling. 

We first state the result without any reference to any operator, recalling 
only the following notation: 

Let n be a positive integer and let pi, . . . , p„ be coprime positive integers. 
Let p = (pi, . . . ,p„). Again let 

G = {CeC*; C''' = 1 for some i} 

and for each C ^ G introduce the index set = {i; C'^' = 1}, whose cardi- 
nality is + 1. Finally recall that m{Q = gcdjpj; i G i^}. 

For N & N and a € N", we define the convex polytope (actually a 
simplex of dimension n — 1) 

V{a, N) = (M+)" n {(xi, . . . , xn); {x, p) = N - (a, p)}. 

Remark 5.6 Notice that V{a,N) is neither integral or Delzant (in the 
terminology of However if p is the least common multiple of pi, . . . , p^ 

then V{pa,pN) has vertices with integral coordinates, but it still not Delzant 
in general. This reflects the fact that V{a,N) is a moment polytope for a 
symplectic orbifold which — except for pi = • • • = p„ = 1 — is not a 
manifold. A 

When a = (ai, . . . , a„) G N" we use the notation al = (ai)! • • • (on)!- 
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Theorem 5.7 For any a G N" we have the asymptotic formula 

(40) 

Moreover 



ao(a, C) = 



Remark 5.8 As we shall see below, the 00(0, C)'s are actually "Weyl 
terms", in the sense that they are the result of some integrals over dif- 
ferent faces of the polytope (or, equivalently, they are phase space integrals 
for some sub-orbifolds). A 

Remark 5.9 If we choose a = 0, the left-hand-side of is just the number 
of integral points of the polytope V{0, N), and the formula becomes: 

00 

#(P(0, N) n Z") = ^ (^~Nj^n(0 N~%{C) + ^(A^~°°) 
CeG 1=0 

where 



-N 



It is known that each factor of C~ is actually a polynomial in N given by a 
Riemann-Roch type formula. See ^] for details on this issue. The leading 
term, obtained with = 1 is 



N 



n-l 



(pi •••p„)(n- 1)! 

The other coefficients can also be obtained as the coefficients of the gener- 
ating function g{X) = U7=ii'^ - ^^0"^ see 1|5|. A 

Proof of the theorem. Let us work in the Bargmann representation. The 
harmonic oscillator we consider is 



j=l - J 
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If a and f3 are multi- indices in N", the differential operator z°'{h-^)^ com- 
mutes with H2{fi) if and only if (a — /3, p) = 0. Thus let us consider the 
symmetric differential operators 

K^ih) := z-{h^r. 

oz 

We shall compute the trace of the restriction of to the eigenspace J^E,h 
in two different ways. The first way is just to do an explicit computation 
in a basis of J^E,h- The second way is to use remark ^.41 in order apply 



theorem 5.5 with g = ld and E = 1. 



1. — A basis of J^E,h (or, more exactly, of the space 'y^jfh ™ 
Bargmann representation) is given by the monomials 

{z^- 7^0and (7,p) =iV} 

where the integer N is defined by the equation E = h{N + |p| /2), and the 
inequality 7^0 stands for 7j ^ 0,Vj. It is straightforward to check that 

, f^'°'r^^^ if 7 ^ a 



otherwise. 



Hence 



7^ q: 
(7.p> = iV 



2. — Formula (|^) is now a simple transcription of theorem 5.5, using 



E = 1 and h = 1/{N + |p| /2), and rewriting the coefficients in order to 
transform the expansion in powers of (|) into an expansion in powers of 
N~'^. This of course does not change the formulas for the leading coefficients. 

Let us compute these leading coefficients. For this we use some homo- 
geneity property in the variable E, so it is best not to set E' = 1 for the 
moment. 

We introduce the standard action-angle coordinates for the harmonic 
oscillator in R^". Let Ij = ^{xj + Cf) = l^jl'^^ and let the angles 6j be 
defined by 




= 

Since d^j A dxj = dij A d9j and the singular set of these action-angle coordi- 
nates is of codimension 2, we can use the corresponding symplectic measure 
on X (M/27rZ)'^ : A • • • A d/n A ^6*1 A • • • A d0n\ as a replacement 

for the symplectic measure on M?^. Let \dE\ denote the pull-back by 
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H2 = Ylj Pj^j of the Lebesgue measure on M. Since the flow of H2 is 27r- 
periodic, the Liouville measure fiE on {H2 = E} is by definition the quotient 
of// by 2TT\dE\ : 

/i = 2TTfiE ® \dE\ . 

Tlie principal symbol of Ka{h) is k^^a = z'^z'^ = H ' fo'^ which we use the 
notation I". Let Sa{E) = f^H^^E}^°'^^E■ 

Since di? A (i/2 A • • • A (i/„, = pid/i A • • • A d/„ we see that 

IJ,E = (27rpi)^"^ \dl2 A • • • A A ^6*1 A • • • A d9n\ 

and 

= {27rr~'pY' [ |d/2 ■■■dln\= E\^\+^-' Sa{l) . 

JpiIl + ---p„In=E 

We conclude by adapting the standard trick used to calculate the surface of 
the unit sphere, namely let 

/•oo 

J = e-^SaiE)\dE\. 
Jo 

J = Sail) / e-^El^l+'^-i \dE\ = Sa{l)T{\a\ + n); 
Jo 

and 



On the one hand 



on the other hand 



e^Pi'i \dh\]---{ I e-P"'" \dln 



00 \ / />oo 



(41) 



Hence 



5 ^^ _ r(ai + l)...r(a„ + l) 



pr^'---pn"+'r(|a|+n)' 
This gives the result for ao{a, 1). For a general ^ the calculation of the 
integral over Mf^ is exactly the same as the integral 5*0, (1) provided we keep 
only indices i G , as follows from the definition of Y"^ in ( |39| ) . This finishes 
the proof. □ 

Remark 5.10 Instead of Ka{h), we could have considered the more general 
operators commuting with H2{K)'- z°'{fi-^)^ , with (a — /?,p) = 0. But using 
the basis z'^ as before it is easy to see that the trace of such operators always 
vanishes as soon as q / /?. As a consequence, this shows that for any C £ G, 

/ z^z'^/i^ = provided (a - /3, p) = 0, a 7^ /?. 

A 
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